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ABSTRACT 



The object of this study was to determine the feasi- 
bility of automatically controlling propeller pitch to 
reduce RPM fluctuations caused by variations in propeller 
speed of advance. The investigation was made ignoring the 
reduced frequency effect and power plant governing. 

A Victory Ship (AP2) was chosen as a prototype and a 
feedback system controlling the pitch of a controllable 
pitch propeller was developed. A non-linear analogue 
computer program approximating the system and the propeller 
characteristics is presented. Electronic difficulties 
prevented obtaining total system response from this program. 
The computer study indicated that a linear analysis would 
be useful. 

The linear analysis indicated that, for the prototype, 
in order to achieve a reduction in RPM fluctuation of 10# 
of that present in the fixed pitch case, a hydralic pump 
consuming over 300 HP would be needed. Systems utilizing 
RPM, shaft torque, and thrust are analyzed. Thrust sensing 
provided the poorest response and in fact increased the RPM 
fluctuation for lower frequencies. 

RPM governing by automatic control of propeller pitch 
was found to be inadvisable unless major design changes 
are made to reduce the spindle torque of presently available 
controllable pitch propellers. 
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INTRODUCTION 

The propulsion systems of ships operating in a seaway 
undergo speed fluctuations caused by variations in propeller 
load. The shaft torque produced by the propeller is a 
function of the flow velocity across the propeller disc 
area. As the ship rolls* pitches, and encounters waves 
this flow velocity fluctuates causing a fluctuation in 
shaft torque which the propulsion system must overcome. 

The propulsion system responds to the fluctuating load by 
allowing fluctuations in RPM. The amount of RPM fluctuation 
present depends on the system Inertia, damping, and the 
effectiveness of the governing, if provided, of the power 
plant. High speed turbine installations with large 
reduction gears experience less RPM fluctuation difficulties 
than do direct drive diesels.* 

The largest variation of flow velocity averaged over 
the propeller disc area, henceforth called speed of advance, 
will probably occur at resonant pitching. Model tests 



* Evans, Harvey J. Conversations before seminar group 
of Department of Naval Architecture and Marine Engineer- 
ing, M.I.T. 24 April 1962. 



indicate that at this condition the variation of speed o$ 
advance could be as high as 0.43 of the steady state 
value. For a geared turbine ship such as the Victory 
Ship (AP2) with a fixed pitch propeller a fluctuation of 
about 6 RPM will result at the frequency of resonant 
pitching, which is about 0.03 cycles per second. Calcu- 
lations leading to these figures are given in Appendix E. 
For systems with less effective Inertia and damping the 
fluctuations would be correspondingly higher. 

Ship's propellers are designed to operate efficiently 
at a designed speed of advance. Their efficiency is 
reduced when Subjected to a constantly varying speed of 
advance. The propulsion system could deliver full power 
when the ship was operating in a seaway if the load could 
be made sufficiently constant. 

A controllable pitch propeller has the ability to 
alter the propeller pitch at any time. If the pitch were 
altered to counteract the effect of the changing speed of 
advance, a constant torque load would be provided to the 
power plant. The efficiency of the power plant and 
propeller might be increased and the entire system could 



Burtls, Evenson M. and J.D. O'Connell. Measurement 
of Ship Model Wake Variation in Rough Water . Naval 
Engineers Thesis, Department of Naval Architecture and 
Marine Engineering, M.I.T. Cambridge, 1959. 






be operated in heavy seas at high power. A literature 
survey has revealed that no reported investigation of 
this possibility has been made. Extensive use of one 
governor systems controlling RPM by altering both pitch 
and throttle setting has been made by the aircraft 
industry in conjunction with turboprop aircraft. 

The purpose of this paper was to investigate the 
possibility of governing propulsion system RPM by auto- 
matic control of propeller pitch. Systems sensing RPM, 
torque, and propeller thrust were investigated. System 
speed governing by controlling the load may or may not 
be more advantageous than conventional governing of the 
power plant. The regulation of some variable other than 
speed might also provide a more optimum overall operation. 
This investigation, however, examines the feasibility of 
speed governing through propeller load control only. 
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II 



GmEBAL PROCEDURE 

Several major analytical assumptions were made at 
the outset. Unsteady flow was tinted using the quasi- 
steady approximation and the reduced frequency effects 
were neglected. The speed of advance was assumed to be 
a steady state value with a sinusoidal fluctuation super- 
imposed thereon. The control systems investigated 
contained a simple one loop feedback with no special 
compensating networks. The power plant was assumed to 
deliver constant torque and no model of the governor was 
developed. 

The standard method of presenting propeller character- 
istics envolves a plot of non-dimensional thrust and torque 
cpefficients as functions of advance ratio, J, and pitch 

ratio P . These plots are very nearly families of straight 

IT 

lines. Por the non-linear study the propeller curves 
were assumed to be families of straight lines. The slopes 
and values of the actual and the assumed curves are the 
same at the operating point. The accuracy of the approxi- 
mation decreases as conditions depart from the operating 
point. The approximations lead to expressions (6) and (13) 



for propeller torque and thrust respectively as derived 
In Appendix A. 

Qp = Kppoo 2 - K 2 V Q 6d - K 3 U> 2 . (6) 

T = K 4 p^ 2 - K 5 V 0 (aJ - Kg CaJ 2 . (13) 

The K's are constants for a given propeller and 
depend on both the propeller curves and the dimensions 
of -the propeller. 

For the linear analysis, the propeller torque and 
thrust were assumed to be linear functions of speed of 
advance, pitch, and RPM. This is the customary incre- 
mental linear analysis and is valid only for small 
variations of the variables about the operating point. 

The equations for incremental variations in propeller 
torque and thrust are derived in Appendix E and are 

Q P = K l* v o + + kV>' (81) 

and 

T = K*Vq + K*/p'+ K*oo’ . (100) 

*. 

The K s are found by assuming a small change in 
the corresponding variable ( + 0.1), and while keeping 
all other independent variables constant, determining 
from the propeller curves the corresponding change in 
the dependent variable. The K* then becomes the 
dependent variable change divided by 0.2 The Troost 
Standard Series for Single Screw Merchant Vessels (B4-55) 
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propeller curves were used throughout. 



No test data providing blade spindle torque directly 
is available. The torque about the blade spindle axis 
was assumed to be composed of hydrodynamic, sliding 
friction, viscous friction, inertia, and centrifugal 
force components. As with the propeller torque and 
thrust coefficients, the hydrodynamic component of blade 
torque can be represented by a family of blade torque 
coefficient curves, K B , where advance ratio and pitch 
ratio are the independent variables. These curves were 
also assumed to be families of straight lines and the 
resulting expression for hydrodynamic blade torque as 
derived in Appendix A is 



forces for a controllable pitch propeller installation 
on a tug,(YTB-502), while underway at various speeds. 

This data also includes details of the pitch setting 
mechanism, ship's speed, centrifugal force components as 
functions of RPM and pitch, and other useful information. 
An approximate wake fraction of 0.164 was assumed based 
on the block coefficient of the tug. Since the tests were 



* Lewis, P.M. "Propeller Coefficients and the Powering 
of Ships," Society of Naval Architects and Marine 
Engineers Paper , Nov. 1951. 



Qbh = KypW ^ - KgV Q u) - K^co^ 



( 21 ) 



Reference (11) contains 




control rod 
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made under steady conditions, the only forces assumed 
acting were hydrodynamic and centrifugal. After 
subtracting the centrifugal component, the control 
rod force was converted to a spindle torque, non- 
dimenslonallzed, and finally plotted as blade torque 
coefficient versus advance ratio and pitch ratio. The 
results are shown in Figure B-3 below and the calcu- 
lations leading to this plot are given in Appendix B. 




FIGURE B-3 

Derived K-g Curves for YTB-502 

- 1 - 



As before stated, the Kg curves were assumed linear. 
The hydrodynamic spindle torque is a strong function of 
the location of the spindle axis relative to the hydro- 
dynamic center of pressure. It was assumed that the 
effect of moving the spindle axis would be to shift the 
K b values up or down the family of K fi curves. Once the 
operating point of the prototype assumed for this study 
was chosen, the operating point of the Kg curves given 
in Figure B-3 was assumed moved to the prototypes 
simply by shifting the K fi values corresponding to the 
tugs operating point to the approximate position of the 
prototype operating point. This can be justified only 
in that such a propeller could be made. The slope and 
spacing of the Kg curves was assumed to remain the same. 

The centrifugal force component of blade torque as 
derived by Boswell and others is 



This expression would be difficult to compute on an 



* Robert J. Boswell, A Method of Calculating the Spindle 
Torque of a Controllable-Pitch Propeller at Design 
Conditions . David Taylor Model Basis Report No. 1529, 
Aug. 1961',' pp. 31-34. 




( 28 ) 
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analogue computer and a simplification given in Appendix 
A leads to expression (35) where the constants K^q and 
are functions of propeller geometry alone. 



The sliding friction component of spindle torque was 
obtained from altering the usual expression for sliding 
friction in a Journal bearing to allow for a moment 
tending to rotate the spindle shaft about a line perpen- 
dicular to the shaft centerline. The resulting expression 



The double sign in expression (38) indicates that the 
direction of the torque must be taken either (+) or (-) 
depending on whether the direction of pitch change is 
toward lower or higher ahead pitch respectively. 

The remaining components of blade spindle torque are 
obtained by direct application of their respective 
definitions. The final expression for blade torque is 



QbC = “ K 10 00 2 ” k h^ 2 p 



(35) 



hydrodynamic 



Q b = Kypw 2 - KQY 0 iO - K 9 u) 



sliding friction viscous friction 




b 



inertia 



centrifugal 



There are numerous variations of hub mechanisms in 



use. All are required to convert motion in the direction 
of the propeller shaft axis into rotation of the blades 
about their spindle axes. As shown in Appendix A, each 
of the four mechanisms most frequently used result in 
approximately the same amount of blade rotation for the 
same input of linear motion along the propeller axis up 
to a pitch ratio of 1.5. The same is true for the force- 
to-torque transmission. For this reason the simplest 
mathematical relation of the four was used throughout. 

The results are approximately applicable to all the 
mechanisms. The expressions used are as follows: 



There are two general types of controllable pitch 
propeller mechanisms. The first, designated herein as 
Design I, utilizes a control rod extending through the 
propeller shaft and is driven by a hydralic piston with- 
in the ship. The second, designated herein as Design II, 
is driven by a power piston housed within the propeller 
hub. The power piston of Design II receives oil from 
a pump within the ship through oil passages in the 
propeller shaft. Design II, as presently available, 
includes a safety spring within the hub that makes the 




(39) 



and 



f CQ = | 



(41) 
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propeller "Fail-safe" by forcing the blades into ahead 
full pitch in the event of a failure in hydralic oil 
pressure. 

The equation of motion of both these designs is 
obtained by equating the oil flow from the pump with 
the sum of the flows into the' mechanism which include 
flow from piston motion, flow from leakage, and flow 
due to compression of the fluid and expansion of piping. 
The final equation of both systems are of the same form 
with the exception of additional terms in the Design II 
expression that are generated by a safety spring force. 
Design I was used for all further study. Schematics of 
both designs are given as Figures A-7 and A-8 in Appendix 
A. . The equation of motion of Design I as derived in 

■1 

Appendix A is 

Input Ram Flow Leakage flow 

Kpx = p 12RA R + ££ 

D Ar 

- k 9^ 2 + (f R + — ^ N ) +_£ cM 

D * tt d R 




+ £_ ( j 1BN + M^R) - Kt nfr) 2 - Kn Ob 2 p 

TTD R R R 




Compliance Flow 
2 



K 7 P OJ ^ _ KqV q CO 



- *2 



R 



R 



60 

R 



-II- 



- <*R 



+ frTjJ N 



bR 



'1 



2 o 

Kip CO - Kjl 60 P 
R R 



D 



£j£ 



V^r r ) 



+ ^ 
t d 



R 



M C R)« 



(42) 



The Victory Ship (AP2 ) was chosen as a prototype to 
determine values of the various parameters. This is a 6000 HP 
geared turbine merchant ship for which considerable data was 
available. The detailed determination of the various 
parameters present in equation (42) above is carried out in 
Appendix B. The propeller blade was assumed to be an 
elliptical, flat plate for the determination of blade inertia 
and centrifugal force components of blade spindle torque. 

The dimensions of the hub were obtained by scaling up the 
design given by Rupp in Reference 11. 

The three remaining parameters to fix were area of the 
ram, the size of the pump, and the hydralic pressure. Since 
the initial cost and maintenance on systems using more than 
1000 psi would be high, this value was chosen as the 
maximum hydralic pressure. The diaximum diameter of power 
piston that could be accommodated within the propeller hub 
has a radius of 15.75 inches. 

Although Design I has the power piston within the ship, 
this was assumed the largest piston size that would be fitted. 
Calculations given in Appendix B show that with the assumed 
hydralic pressure of 1000 psi and this ram area there is still 
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a likelihood of stalling the mechanism (by opening relief 
valves) even with this large size power piston. The 
maximum combinations of blade centrifugal, hydrodynamic, 
and sliding friction when acting together require 1000 in 2 
of piston area. The assumed piston area is 775 in 2 . 

Larger piston sizes would require larger pump capacities 
which are excessive for this piston area. 

As shown in Appendix E, the time constant of the pitch 
control mechanism is 

^p= fis — . (so) 

KpTTDM 

The effectiveness of the system in regulating RPM 
fluctuations requires that this time constant be small. 

The larger pump capacities correspond to larger Kp. Assuming 
that the marine industry would not accept a device consuming 
more than 5 % of the installed main plant horsepower, a 
maximum pump size of 300 HP was assumed. Kp was then 
determined from manufacturers data. 

The viscous friction forces were obtained from an 
analysis of transient spin test data taken by Rupp. This 
data consists of a plot of control rod forces versus pitch 
for the operating conditions of going from full ahead to 
full astern and vice-versa. This data is reproduced as 
Figure B-7 below. The thrust forces causing high values 
of sliding friction were assumed zero for this test. The 
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FIGURE B-7 

YTB-502 Control Rod Force for Spin Test 
in Air From Data by Rupp 

long region of both transient curves that are parallel 
to the steady state condition are assumed to correspond 
to the condition of maximum pump output and constant piston 
speed. The deviation from the steady state curve in both 
cases is approximately 3 x 10^ lbs. This control rod 
force was converted to spindle torque and scaled up by the 
ratio of propeller diameter squared since viscous friction 
is proportional to surface area. The spaed of rotation 
was assumed constant and equal to the total angle divided 
by total time. 
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The substitution of the various parameters into 
equation (42) revealed that several, of the terms would 
be very small in comparison to others. The coefficients 
of each independent variable was examined and those 
terms having coefficients of less than 1 % of those re- 
maining for the same variable were eliminated. A 
further simplification was made by setting to zero all 
terms including the first derivative of shaft speed. 

This assumption is based on the knowledge that propeller 
RPM does not vary to a large extent even without regu- 
lation and should vary less with regulation. The 
resulting simplification of equation (46) is 



iux = p 

TT D 



(K 7 -Kh)PUJ 2 _ K 8 v 0 60 



R 



R 





(65) 
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Ill 



PROCEDURE AND RESULTS OF 
THE NON-LINEAR ANALOGUE COMPUTER STUDY 



A. Procedure 

Equations (6), (13), and ( 65 ) given above envolve 
non-linearities in the form of products of Independent 
variables and a discontinuous sliding friction term. If 
saturation of the pump is also allowed, another discon- 
tinuity is introduced which would limit the value of pump 
output, kpX, to some finite value. All of these non- 
linearities can be programed on the analogue computer 

that was available. The program includes four multipliers 
0 p . 

to generate^ , p(*j , A and v Q u) . It also Includes a 
switching mechanism controlled by the derivative of pitch 
to approximate sliding friction. 

The program was arranged so that shaft torque and 
thrust were generated in order to test the final system 
sensing these variables. 

The program is shown in Figure C-2 below. The para- 
meters determined in Appendix B were then inserted into 
the program and voltage scales were assumed for each 
variable. The gains of various loops were then adjusted 
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to provide signals that would not saturate the various 
amplifiers and multipliers of the computer. The details 
of this scaling is given in Appendix C. 

The input to the system was taken to be the voltage 
analogue corresponding to speed of advance, V Q , and 
consisted of a steady positive value with a sinusoidal 
voltage of amplitude &V 0 provided by a signal generator. 
The output consisted of a voltage analogue for shaft 
speed, to , and both input and output were recorded on a 
two channel Sanborn recorder. It was hoped that data over 
a range of frequencies could be obtained in this manner 
for the fixed pitch case as well as the system operating 
in closed loop sensing shaft speed, torque, and thrust. 

B. Results of the Non-Linear Study 

The frequency response of the fixed pitch system was 
obtained. The data is recorded in Appendix C and is 
shown graphically as a plot of the magnitude of the ratio 
of output to input in decilogs in Figure I below. The 
strong stabilizing effect of sliding friction was very 
apparent. When the signal corresponding to sliding 
friction feedback was removed with a sinusoidal AV 0 
input of amplitude 6 ft/sec the second derivative of 
pitch, increased by a factor of 80. The feed forward 
path containing Vo<aj had no effect on p even for very 
large values of &V 0 . 
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FIGURE I 

Fixed PiVc^ Frequency Response 



Values generated by the computer for propeller 
torque, thrust, and blade torque were within 5$ of the 
calculated values. The time constant of the system with 
fixed pitch was observed by injecting a step in engine 
torque and recording the decay in shaft speed variation. 
The time constant was observed to be approximately 3.28 
seconds. The decay to steady state required about 10 
seconds (approximately 3 time constants). 

Electronic difficulties, the details of which are 
given in Appendix D, prohibited obtaining data for the 
control system operating in closed loop. 

In all cases the node containing p was very small 
and consisted primarily of computer generated noise. 
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IV. 

PROCEDURE AND RESULTS OF THE LINEAR STUDY 
A. Procedure 

The results of the non-linear study Indicated that 
valid further simplifications could be made in equation 
(46). Eliminating the terms containing Vto and p as 
suggested by Chapter III above, equation (46) reduces to 



Equation (67) includes ram flow and leakage terms 
only. Since the linear analysis is valid only for small 
variations around the operating point, the leakage was 
assumed constant and is equivalent to an effective 
reduction in pump output Kp//. As shown in Appendix E, 
equation (67) eventually reduces to a formula whose 
frequency domain transfer function is 



where y’(s) = Laplace transform of the incremental 




(67) 




p'(s) = £ZM 

y'(8) TpS +1 



mechanical input to the pitch control mechanism. 

Using the linear approximations for propeller 
response before mentioned, and assuming the power plant- 
propeller system to be a first order system containing 
inertia and damping, the complete systems for sensing 
RPM, torque, and thrust were drawn in block diagram 
form as shown in Figure II. 

The gains G^ and Gq for the shaft speed and shaft 
torque sensing systems respectively were chosen so that 
a damping ratio of 0.4 was assumed by the closed loop 
system. Higher values of loop gain reduce the low 
frequency response but cause higher peaking or ampli- 
fication at a higher frequency. This peaking is apparent 
at about 2 radians per second as shown in Figure E -6. The 
gain for thrust, G^, was taken to be a very low value 
only to demonstrate that any positive value of G^ is not 
beneficial in regulating RPM at the lower frequencies. 

B. Results 

System gain is defined as the ratio of shaft speed 
variation for a given variation in speed of advance, ^ * . 

V 

The system gain for each of the three linear systems 
pictured in Figure II was obtained and plotted by the 
standard methods. 

The results are shown in Figure E-6 where system gain 
is plotted in decilogs. 
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(a) Regulation by Shaft Speed Sensing 




(b) Regulation by Shaft Torque Sensing 




(c) Regulation by Thrust Sensing 



FIGURE II 

Linear Approximation to Complete Systems 
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DISCUSSION OP RESULTS 

As shown in Figure E-6, the fixed pitch response 
determined on the analogue computer was in good agreement 
with the linear analysis for frequencies below 0.10 rad/sec. 
It is of interest that the input levels used to obtain this 
response curve are extremely high. The minimum value for 
the fluctuating component of the speed of advance is three 
volts which corresponds to thirty feet per second. The 
steady state value of speed of advance is only 16.7 ft/sec. 
The maximum value of input reached 40 volts corresponding 
to 400 ft/sec. These high values of input were necessary 
to obtain an output above the noise level. The results 
therefore indicate that the linear analysis should be 
reasonably valid even for moderate fluctuations around the 
operating point. 

Had it been known that the terms containing p and V in 
equation ( 65 ) would be insignificant, an analogue eliminating 
these paths would have been developed. As explained in 
Appendix D, hum in the *p path made the obtaining of data 
for the closed loop system impossible. The simpler analogue 
would not have this low signal level path and might have been 



more successful. Unfortunately time was not available 
to reprogram. 

The simplifications leading to equation (65) as 
shown in Appendix C, call for ignoring the term 

— ( f R + f r T -L^ M ) 

dt bR 

This friction term represents an oil flow and since 
it is discontinuous at p = o, corresponds to an impulse or 
hydralic shock. This hydralic shock would very likely be 
damaging to relief valves and other parts of the system. 

The introduction of an accumulator to cushion this shock 
would introduce additional compliance into the system, but 
as shown in Appendices C and D compliance is relatively 
unimportant. There are no compliance terms remaining in 
the final linear analysis. 

The simplifications of expression (46) given in 
Appendix C show that inertia loading is relatively unimpor- 
tant. Design I has a larger inertia loading because of the 
relatively larger mass of the control rod compared to the 
smaller power piston of Design II. The area of the high 
pressure port delivering full flow to the power piston of 
Design II would have to be 9.1 square inches to match the 
output of the variable displacement pump at the same pressure 
Since inertia loading is relatively unimportant and since 
providing a large oil port area would be a difficult engineer 
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ing problem, there Is little advantage to Design II for 
this application. 

The transfer functions used in comparing the various 
systems are stable for all frequencies. These transfer 
functions are based on numerous assumptions that neglect 
many of the lags and storage elements actually present. 

They also do not consider nonlinearities such as sliding 
friction, strict ion, saturation, and the nonlinearities 
present in the propeller characteristics themselves. 

There is therefore a possibility that the actual hardware 
might exhibit natural frequencies corresponding to 
resonances in the remainder of the propulsion system. 

The lowest frequency of torsional vibration resonance for 
the Victory Ship (AP2), as measured and reported in 
Reference 16, is 28 cycles per second or 175 radians per 
second. 

With the restrictions of a maximum pressure of \000 
psi. hydralic pressure, a maximum pump power of 300 HP, 
and the approximate dimensions of present day controllable 
pitch propellers, the improvement in shaft RPM regulation 
over a fixed pitch propeller installation using automatic 
pitch control corresponds to a reduction in RPM fluctuations 
by a factor of nine. The improvement in propeller and 
power plant efficiency would not be great enough for this 
small an increase in RPM regulation to warrant the 
installation of such a large and complex system. 
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Several improvements in the system analyzed might 
reduce the size of the equipment necessary or improve 
the system response. The raising of the hydralic 
pressure would reduce the size of the piston necessary 
and thereby reduce the capacity of flow necessary from 
the pump. This would improve the system by reducing 
the important pump-piston time lag. Increasing the 
horsepower of the pump also improves the system by 
decreasing the pump-piston time lag. Compensation in the 
form of lead or lag networks in the feedback loop would 
delay the peaking present in the proposed systems and 
allow a higher closed loop gain and hence a higher degree 
of regulation. 

It is unlikely that the above mentioned changes 
would result in improvements^ system performance to any 
large degree. If spindle torque could be greatly reduced, 
the size of the power piston could be reduced and the 
entire system could be designed wi^bh better response 
and be made smaller. As shown by equation (38), spindle 
torque consists of many components, each of them large. 

By increasing hub size the spindle bearing could be made 
larger and the lever arm of the hub mechanism increased. 

The application of non-friction bearings might reduce 
this factor although roller and ball bearings are presently 
considered unacceptable for naval service due to shock 
susceptabllity . If the propeller were designed to operate 
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at one RPM, the hydrodynamic and centrifugal components 
might be reduced by careful design. These significant 
changes might make automatic control of pitch to remove 
fluctuation in RPM feasible. 

As demonstrated by Figure E-6 there is no advantage 
in sensing either torque or thrust in lieu of RPM when 
the variable being regulated is RPM* The improvement 
in regulation offered by RPM sensing is 9.7 decilogs 
which corresponds to a reduction in fluctuation by a 
factor of 9. 
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VI 



CONCLUSIONS 



For geared turbine qr geared diesel propulsion plants 
used by the U.S. maritime industry, automatic control of 
prqpeller pitch to reduce wave excited RPM fluctuations 
would be excessively expensive in horsepower, space, and 
complexity for a small gain in RPM regulation. 

Significant reductions in dynamic blade spindle torque 
through design changes would be necessary to make this 
scheme practical. 



VII 



RECOMMENDATIONS 

Automatic pitch control investigation requires a 
knowledge of spindle torque. Very little is known In 
this area. It is recommended that future effort be con- 
centrated on the determination of dynamic spindle torque. 

A more useful approach to the problem of automatic 
pitch control would include the presence of engine 
governing with a view toward developing a system that, 
with one governor would receive a bridge signal and 
control pitch and engine throttle setting to achieve the 
thrust or ships speed desired by the operator. 



VIII. 

APPENDIX 



APPENDIX A 



Derivation of System Equations 

1. The Power Plant-Propeller System: 

Figure A-1 represents a schematic of the propeller- 
power plant system. 

The equation of motion of this system is 

QE - QP - <£ e w = (Jip + Ji E ) & > (1 ) 

where 

Q E = developed engine torque, lb-ft 
Qp = absorbed propeller torque, lb-ft 

i 

Jip = effective polar mass moment of inertia of the 

propeller (includes entrained water), lb-ft-sec^ 

i 

Ji E = effective polar mass moment of inertia of the 

engine and gearing referred to the propeller shaft, 
lb-ft-sec 2 

= mechanical damping in the power plant and gearing, 
ft-lb-sec. 

Q s = shaft torque, ft-lbs. 

Co = shaft rotational speed, radians / second 

• O 

to = shaft rotational acceleration, rad/sec . 

The torque delivered by the engine is 

Qe • 

The shaft torque expression is 

Qs = Q e - oC-gto - JljA) . (2) 

Curves of Kq for most propellers are nearly linear over a 
wide range of adrance coefficient and pitch ratio. Figure A-2 
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FIGURE A-I 

The Propeller-Power Plant System 




FIGURE 2 

Assumed Linear Torque Coefficient (Kq) Curves 
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represents an approximate general set of Kq curves assuming 
that they are completely linear. 

The equations of these lines are 

(3) 



f = + + hQ 



where 



K C 

P 

D 



Q 

= at 



= propeller pitch, ft, 

= propeller diameter, ft, 



trn q = slope of the Kq curve at the operating point of the 
true Kq curve and assumed constant for this analysis, 
non-dimensional 

J = propeller advance coefficient = ^ 

Vo = speed of advance for the propeller, ft/sec 

n = propeller speed, rev, /sec 
A ± 



— ratio 



, non-dimensional, assumed constant 



. A <.£2 

hQ = / f/D intercept of the Kq=0 line. 

^ = water mass density 1 .024^ 

When nn q , bQ, and hQ correspond to those at an operating 
point taken from actual propeller curves, the response of the 
propeller to changes around the operating point will be 
accurate within the limitations of a steady flow analysis. As 
the propeller operates farther from this point the assumption 
of linearity becomes less valid, but is considered sufficiently 



1 . This unit for mass density for water was taken to be 
consistant with Troost's propeller curves used to 
determine coefficients of equation (3). 
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accurate for this study. 

Solving equation (3) for Qp, 



Qd = 'P ? 2 d 4 - K)qVo P^D 4 - h Q Pn 2 D 5 
b Q b Q b Q 


t 


(4) 


Since n = — , 

2TT 

equation (4) can be written in the form 




(5) 


Q p = K^fu> 2 - K 2 VoOJ- K 3 CD 2 ; 

where 




(6) 


K 1 " 

W^Q ’ 




(7) 


K 2 = m 3f£l. , 

2 TTbQ 




(8) 


and K x = ho P D 5 
4tT 2 bQ 




(9) 



Substitution of equation ( 6 ) into equation (1 ) completes 
the derivation of the propeller-power plant system equation. 

Q e - K if 00 2 + K 2 VoO> + Kj60 2 -cC e O) = ( Jip + Ji E ) OJ (10) 

The thrust coefficient for most propellers can also be 
considered linear with respect to 'P/D and J. Again this is 
accurate for small changes around an operating point and less 
accurate for large departures from this point. Figure A-3 
represents the linear approximation to the actual K T curves. 

The procedure for obtaining thrust as a function of Vo, go, 
and p follows in the same manner as used in obtaining the 
similar expression for Qp. 

From the linearized K T curves 

£ = J + b T K T + h T „ (11) 
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FIGURE A-3 

Linear Approximation to K T Curves 




Linearized Hydrodynamic Blade Torque Coefficient Curves 
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( 12 ) 



Using the definition of K T , 
K T= ^ n^D 4 » 



and solving for T . 

T = K li'fUO 2 - K 5 Vo^) - k 6 uo 2 



where 



K 4 = 

k 5 = 



4 TT £ ^r. 

P D' 






and Kg = 
With b = 



ZTT Jr T 

hr PP* 

A K-r 

^n T = slope of K T curve, 
and h.p = J? intercept of Kq = 0 line. 



(13) 

(14) 

(15) 

( 16 ) 

(17) 



2. Spindle Torque 

The equations of motion of any servo system are a function 
of system loading. Loading in this case consists of blade 
spindle torque. Torque about the spindle axis can be viewed 
as composed of hydrodynamic, sliding friction, viscous friction, 
inertia, and centrifugal force components. 

These components are derived separately. The equation 
for spindle torque is 

Qb = q bh + Qbs + Qbv + Q BJ + Q BC ( l8 ) 

where 

Q b = total spindle torque that must be overcome by the 
pitch changing mechanism, lb-ft 

Qbr = hydrodynamic spindle torque, lb-ft 

Q B g = sliding friction torque, lb-ft 
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Qbv = viscous friction torque, lb-ft 
QBJ ~ Jig ^ b» lb-ft 

Q B q = centrifugal force torque from the mass of the 
propeller metal, lb-ft. 

Very few curves of spindle torque are available. Hydro- 
dynamic spindle torque is a strong function of the spindle 
axis position relative to the center of hydrodymic 
pressure. Tests made at M.I.T. show that curves of K B , where 
Kg is defined as 



are very nearly linear functions of advance coefficient and 
pitch ratio just as is As with Qp linear curves for Kg 

are assumed as shown in Figure A-4, and the expression for 
Qbh is derived based on this assumption. 

For these curves a spindle torque that tends to increase 
pitch is considered a positive torque. 

Proceeding as before. 




(20) 



and 



Q BH = K 7 f>u>2 - KgVo - K 9 o> 2 



(21) 



where 




( 22 ) 




Conversation with S. Curtis Powell, March 2, 1962 
M.I.T. 



- 39 - 



( 23 ) 



Kq = P P 4 
2 TT b fi 



and Kg = h B P 1)5 
4 tt^ b 



(24) 



The spindle Is not a journal loaded In the normal 
manner. Propeller thrust acts at a considerable moment 
arm tending to rotate the splnde about an axis perpen- 
dicular to Its centerline as shown In Figure A-5. 

If the bearing were normally loaded the friction torque 
would be 

fPr , lb-ft j 

where 

f = coefficient of journal friction 
P = bearing load, lb 
r = spindle radius, ft. 

Tj_ = total thrust on blade perpendicular to the spindle 
axis, lb. 

An approximate expression for spindle sliding friction is 



Ji = lever arm from center of thrust on the blade to the 
center of the spindle bearing, ft. 

b = spindle bearing length, ft. 

N = number of blades 

Sliding and static friction are assumed to be the same. 
Blade viscous torque can be expressed as 



Qbs = fr 

b 



(25) 



where 



QbV = «£ B CaJ) B Kl 



( 26 ) 
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Spir»de bearings 



FIGURE A-5 

Schematic Indicating Blade Sliding Friction (Qg S ) 




The Relationship Between Pitch and Blade Angle 
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